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ABSTRACT 

It is shown that Fermi acceleration at an ultra-relativistic shock wave cannot operate on a particle 
for more than 1 1/2 Fermi cycle (i.e., u ^ d ^ u — > d) if the particle Larmor radius is much smaller 
than the coherence length of the magnetic field on both sides of the shock, as is usually assumed. 
This conclusion is shown to be in excellent agreement with recent numerical simulations. We thus 
argue that efficient Fermi acceleration at ultra-relativistic shock waves requires significant non-linear 
processing of the far upstream magnetic field with strong amplification of the small scale magnetic 
power. The streaming or transverse Weibel instabilities are likely to play a key role in this respect. 
Subject headings: shock waves - acceleration of particles - cosmic rays 



1. INTRODUCTION 

Fermi acceleration of charged particles bouncing back 
and forth a coUisionless shock wave is at the heart of a 
variety of phenomena in high energy astrophysics. Ac- 
cording to standard lore, this includes the acceleration 
of electrons at gamma-ray bursts internal/external rel- 
ativistic shock waves, whose synchrotron light is inter- 
preted as the prompt/afterglow radiation. 

However the inner workings of Fermi acceleration at 
relativistic shock waves remain the subject of intense 
study and debate, even in the test particle limit. It 
has been argued that a universal energy spectral index 
s ~ 2.2 — 2.3 should be expected (e.g., Bednarz & Os- 
trowski 1998, Achterberg et al. 2001, Lemoine & Pel- 
letier 2003, Elhson & Double 2004, Keshet & Waxman 
2005), which would agree nicely with the index that is 
inferred from gamma-ray bursts observations. Yet re- 
cent numerical simulations that include more realistic 
shock crossing conditions have indicated otherwise (e.g., 
Lemoine & Revenu 2006, Nicmicc & Ostrowski 2006), 
so that the situation is presently rather confuse. Un- 
fortunately, numerical Monte-Carlo simulations of Fermi 
acceleration, although they remain a powerful tool, do 
not shed light on the physical mechanisms at work. 

In the present Letter, we offer a new analytical discus- 
sion of Fermi acceleration in the ultra-relativistic regime. 
We rely on the observation that, under standard assump- 
tions, the Larmor radius of freshly injected particles 
is much smaller than the coherence length Zcoh of the 
upstream magnetic field, and we integrate the equations 
of motion to first order in the quantity rj^/lcoh (Section 
2). We thus find that a particle cannot execute more 
than ll/2u^d^u cycles through the shock be- 
fore escaping downstream, hence Fermi acceleration is 
mostly inoperative. This result is related to the very 
short return timescale in the ultra-relativistic regime: 
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over its trajectory, the particle experiences a nearly co- 
herent and mostly transverse magnetic field, hence the 
process is akin to superluminal acceleration in regular 
magnetic fields discussed by Begelman & Kirk (1990). 
We also show that our analytical predictions agree very 
well with numerical simulations. Finally, we argue that 
the (expected) non- linear processing of the magnetic field 
at ultra-relativistic shock waves and, in particular, the 
strong amplification of small-scale power, may be the 
agent of efficient Fermi acceleration (Section 3); we sug- 
gest new avenues of research in this direction. 

2. ANALYTICAL TRAJECTORIES 

2.1. Field line curvature 
The upstream magnetic field consists of a regular com- 
ponent Bq and a turbulent component SB: B = Bq + SB, 



with (B^) = B'^ + {SB'^). The turbulence is defined in 
the wavenumber range fcmin < k < fcmax by its power 
spectrum S{k) oc fc~", which is normalized according to: 
J d'^kS{k) — {SB'^). Downstream and upstream mag- 
netic fields are related to each other by the MHD shock 
jump conditions (see e.g.. Kirk & Duffy 1999). In the 
case of 7— ray bursts external shocks, the inferred frac- 
tion of energy density stored in magnetic turbulence is 
of the order of a percent (Gruzinov & Waxman 1999) so 
that the magnetic field can be considered as passive. In 
this limit, one finds: — , Bj_\^ = i?sh^j.|u; 

= B ■ z is the component of the magnetic along 

the shock normal z, B± is the projection of B on the 
shock front plane {x,y), and subscripts |d (resp. |u) 
indicate that the quantity is measured in the down- 
stream (resp. upstream) plasma rest frame. The quan- 
tity Rsh = rsh|u/3sh|u/(rsh|d/3sh|d) is the proper shock 
compression ratio, expressed in terms of l3sh\u (resp. 
Psh\d) the shock velocity measured in the upstream (resp. 
downstream) rest frame and the corresponding Lorentz 
factor Fghiu (resp. Fghid). In the ultra-relativistic limit 

(Fshiu > 1): ^sh ^ r,h|u%/8 > 1. 

Hence, to an error ~ 0{l/rsYi\u) on the direction of B, 
it is a good approximation to consider that the magnetic 
field lies in the transverse {x, y) plane downstream of the 
shock. The magnetic field at a given point Tq on the 
shock surface can be written in both downstream and 
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upstream reference frames as follows: 

B\d{ro) ~ B±\d cos{(I)b)x + B±\d sm{(pB)y , 

B|u(ro) = Sx|u cos{(f>B)x + B_L|u sm{<j)B)y + (1) 

The phase 4>b is invariant under the Lorentz transfor- 
mation from downstream to upstream; this observation 
plays a key role in the discussion that follows. 

We assume for the time being that the Larmor radius 
rL of the test particle is much smaller than the coherence 
length of the magnetic field Icoh, if a > 3, Zcoh ~ 1/fcmin 
as the magnetic power is distributed on the largest spa- 
tial scales. Then, since the typical u ^ d ^ u cycle 
time through the shock is of order O (rL/Fghiu) 'coh 
(Achterberg et al. 2001, Lemoine & Pelletier 2003, 
Lemoine & Revenu 2006), in a first approximation one 
can neglect the magnetic field line curvature over the tra- 
jectory of the particle. 

This approximation may be justified as follows. Con- 
sider a particle moving over a length scale I ~ ''L/rghlu ^ 
/coh- The radius of curvature TZyi of the magnetic field 
on scales larger than / can be calculated as: 



velocity vector as: 



'^>i — 



{B-V)B 



1/2 



(2) 



Assuming that a > 3 (or Zcoh ~ 1/fcmin) and decom- 
posing 5B in a Fourier series, one finds that indeed, 
^/7^>^ - {SB/B)~^ (V'coh)/""^^/^ < 1- Hence the large 
scale component is approximately uniform over a length 
scale It is safe to neglect the magnetic power on scales 
smaller than I since SB'^^, - SB'^{l/lcoh)°'~^ < SB'^, the 
latter being comparable to the large scale component. If 
a < 3, similar conclusions apply, since the assumption 
Th ^coh translates into tl ^ ^mixi which means that 
the particle only experiences a smooth large scale mag- 
netic field. 

Now, if the magnetic field is approximately regular over 
the path of the particle in a u d — > u cycle, Fermi ac- 
celeration in the ultra-relativistic regime becomes similar 
to superluminal acceleration in a fully regular magnetic 
field, which is known to be inefficient (Begelman & Kirk 
1990). In the following, we extend the discussion of these 
authors and compare the predictions to numerical simu- 
lations of particle propagation in realistic turbulence. 

2.2. Analytical trajectories 
Upstream. The equation of motion reads: 



d^ 
dt 



PxB 
B 



(3) 



with P the velocity of the particle and JIl = c/rL the 
Larmor frequency. Shock crossing from downstream to- 
ward upstream requires > Psh\u with p'^^"^^ the ingress 
component of the velocity along the shock normal. Hence 
C'(l/Fsh|u)- By working to first order in l/F^hiu, 
Achterberg et al. (2001) were able to obtain analytically 
the particle trajectory and its direction at shock recross- 
ing u — > d. Assuming that (j)B = so that the transverse 
component of B lies along x, one obtains the outgoing 
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Downstream. There we must proceed differently as the 
return timescale ^ 0{rh/c) can no longer be treated as 
a small quantity (rL is evaluated in the downstream rest 
frame; Lemoine & Revenu 2006). 

We may assume that (j)B — since the phase is pre- 
served by the Lorentz transformation; furthermore shock 
compression results in a magnetic field essentially ori- 
ented transversally to the shock normal. To this order 
of approximation, the trajectory along the shock normal 
reads: 

ni^z{t) = sin((/.i) [cos(riLi) - 1] + P^, sin(f7Lt) , 

(5) 

where (f)i denotes the phase of the ingress velocity vec- 
tor in the {x,y) plane: = p'f\ cos{(j)i) , 



/J^j^"*; sin(0i). The shock front follows the trajectory: 

•2^sh|d(^) — Psh\dt and return to the shock will occur if 
and when: 



h) = 9{i) 



/3sh|d^- /3||,iSin(£) 



P^f\ [cos(t) - 1] 



(6) 



with i = Cli,t. The function g{i) diverges toward — oo for 
t ^ 0, 27r and its derivative is monotonous in the interval 
t g]0, 27r[. Hence return to the shock can occur if and 
only if the maximum of g{t) exceeds the value sin(0i). 
Note also that g{t) is always negative since by assumption 
P^l^l ^ /3sh|d- Therefore a necessary condition for return 

to the shock front is <j)i G [— tt, 0], or equivalently P^'^] < 0. 

Once the time t^ of shock return has been determined 
(numerically), the outgoing velocity vector can be de- 
rived from the solutions to the equations of motion: 



Pf, cos(17Lid) 

pf} cos{n^td) 



- Pi"^] sin(f^Lid) 
/3j^] sin(f^Lid) . 



(7) 



2.3. Mappings: downstream to upstream and vice-versa 

Equations Q and ((Tj) define mappings from the ingress 
to the egress angles on either side of the shock. The 
ingress angles in one rest frame are related to the egress 
angles in the other rest frame by the Lorentz transfor- 
mations: 



a(") 



Pi 



(u) 

II, i 



1 -I- Pii^lProl 



4") 



(d) 



^11, f 



1 



(8) 



where P^ci = {Psh\u - /3sh|d)/(l - /3sh|u/3sh|d) is the rela- 
tive velocity between the upstream and downstream rest 
frames. Using these mappings and transformations, one 
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can follow the trajectory of a particle. Since the cycle 
time is of order ri^/T sh ^ Icoh, it is reasonable to assume 
that (j)B = remains constant from one Fermi cycle to 
the next. 

Now, in Section 2.2.2, we argued that /Sy'^l < is a 
necessary condition for the particle to be able to return 
to the shock. However, as the particle travels upstream 
and exits back toward downstream, its outgoing velocity 

is given by Eq. Q and it can be shown that l3^\ > 
irrespectively of the upstream ingress angle. In effect, 
for a given /3^"? , the quantity /3^"| is minimal when /3^\ 
is maximal, i.e. when = 1 — Psh\u ~ f^y^l^ ■ Then 

the final = (-^"i + ^If^lj^ll) /2 > 0. The minimum 

is then 0, which corresponds to a particle entering up- 
stream along X (tangentially to the shock surface), i.e. 

f^z^l ~ Ah, = 0. Hence, if a particle that travels 
downstream is able once to return to the shock, it will 
not do so in the subsequent cycle. 

A quantitative assessment of this discussion is shown 
in Fig. n which presents the locii of ingress and egress 
velocity vectors in the {x, y) plane as seen in the up- 
stream rest frame. The blue area shows the region of 

egress and /3^"f (equivalently ingress as seen from 
downstream) for which the particle is bound to return 
to the shock. The green circles show the ingress 
and /9^") of a particle that crosses toward upstream. The 

various circles correspond to different values of upon 
entry; the radii of these circles are bounded by the shock 
crossing condition /3^"; > /3sh|u- Finally, the red kidney 
shaped forms map these ingress upstream velocities into 
the egress velocities, according to Eq. Q). The fact that 
these kidney shaped forms do not overlap anywhere with 
the blue area confirms that at most one and a half cycle 
u^d^u^dis permitted. 

2.4. Comparison with numerical work 

The previous discussion relies on several approxima- 
tions, most notably that the field lines can be considered 
as straight over the trajectory of the particle. Compari- 
son of the previous results with numerical simulations of 
particle propagation in refined descriptions of the mag- 
netic field are best suited to assess the error that results 
from these approximations. Figure |21 which shows the 
contour plot of the return probability defined as a func- 
tion of Px. i and i, can be directly compared to the blue 
area of Fig. ^ Indeed, the agreement is excellent. The 
parameters of the simulations whose results are shown 
in Fig. 121 are as follows: rL/Lmax = 7 x 10~^, a = 11/3 
(Kolmogorov turbulence), Bq — Q (pure turbulence) and 
Tshiu = 38. The numerical procedure used to follow the 
particle trajectory has been described in Lemoine & Pel- 
letier (2003) and Lemoine & Revenu (2006). 

This discussion also explains the results of recent 
Monte-Carlo simulations. For instance, Niemiec & Os- 
trowski (2006) report that Fermi acceleration is ineffi- 
cient in the ultra-relativistic regime for upstream Kol- 
mogorov turbulence; their simulations indicate very steep 
spectra if any, in good agreement with the present discus- 
sion. In contrast, other studies of Fermi acceleration ob- 
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Fig. 1. — Mapping from downstream to upstream and back to 
upstream as measured in the upstream rest frame, in the plane 
transverse to the shock front; Px and (3y are the velocity compo- 
nents in this plane (note the enhancement by on each axis). 
The solid blue area shows the region of egress upstream coordinates 
which permits the particle to return to the shock from downstream. 
The green circled area shows how the original downstream parti- 
cle population maps upon entering upstream, and the red kidney- 
shaped region shows the mapping of this population on exit from 
upstream. 
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Fig. 2. — Same as Fig. [Tl but using results from the numeri- 
cal simulation of the trajectories of ~ 10^ particles in shock com- 
pressed turbulence (see text for details). Each contour represents a 
drop in return probability by a factor 2.4; the lighter the shading, 
the higher the return probability, black corresponding to zero. The 
yellow circles and red kidney-shaped curves are the same as those 
shown in Fig. [Tl 



tain powerlaw spectra of various spectral indices. How- 
ever, one can check that these latter studies have, one 
way or another, either assumed an isotropic downstream 
turbulence, or implicitly marginalized over the angle be- 
tween the particle trajectory and _Bx at shock crossing, 
which amounts to picking (pB at random in each half- 
cycle. In the light of the above discussion, it is then easy 
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to understand why Fermi acceleration seemed efficient in 
these studies. 

Strictly speaking, our results do not apply to scale in- 
variant turbulence, i.e. a — 3. However, results of nu- 
merical simulations for this particular case are similar to 
those shown in Fig. |21 the return probability is non-zero 
everywhere but it is ten times lower in the kidney-shaped 
region than in the negative j3y^i area. This suggests that 
quite steep powerlaw spectra should emerge from Fermi 
acceleration in such turbulence. 

3. DISCUSSION 

Fermi acceleration is thus inefficient at ultra- 
relativistic shock waves if the Larmor radius tl of in- 
jected particles is much smaller than the coherence length 
'coh of the turbulent magnetic field on both sides of the 
shock. This does not mean that Fermi acceleration is 
bound to fail. In particular, if tl ^ Icoh-, powerlaw spec- 
tra must emerge as the memory of the magnetic field 
direction (in the transverse plane) at shock crossing is 
erased during propagation in small scale turbulence (as 
we have checked numerically). The final value of the 
spectral index will depend on the transport properties of 
the particle in this small-scale turbulence. It is difficult 
to probe this regime using numerical simulations as inte- 
gration timescales become large (the scattering timescale 
3> r-^/c). On analytical grounds, one expects s ~ 2.3 if 
the small scale turbulence is isotropic downstream of the 
shock wave (Keshet & Waxman 2005); if it is anisotropic, 
one ought to expect a different value however. 

Although Fermi acceleration could operate efficiently 
on high energy initial seed particles with tl ^ ^coh, the 
abundance of such particles is generally so low in real- 
istic astrophysical shock wave environments that the in- 
jection efficiency would be extremely small (see Gallant 
& Achterberg 1999 for instance). 

In the above context, the interpretation of the after- 
glow emission of 7— ray bursts as the synchrotron radia- 
tion of electrons accelerated at the ultra-relativistic ex- 
ternal shocks (Fshju ^ 100) becomes particularly enlight- 
ening. The success of this model indeed requires both 
efficient Fermi acceleration as well as very significant 
amplification of the interstellar magnetic field (Gruzi- 
nov & Waxman 1999). Hence it is tempting to tie these 
two facts together and to wonder whether this non-linear 
MHD processing could not be the agent of efficient Fermi 



acceleration. 

One proposal discussed so far is the transverse Weibel 
instability which could possibly produce sufficiently 
strong magnetic fields on very small spatial scales 
- 105cm(r,h|u/10)"^/^(ne/lcm-3)-i/2 (Medvedev & 
Loeb 1999); there is however ongoing debate on the life- 
time and strength at saturation of the magnetic field 
(e.g., Wiersma & Achterberg 2004; Lyubarsky & Eichler 
2006). Nonetheless, such a small-scale turbulence should 
result in powerlaw spectra of accelerated particles; albeit 
the value of the resulting spectral index is not known. 

Recently, it has been suggested that the generaliza- 
tion of the streaming instability to the relativistic regime 
could amplify the magnetic field to the values required 
by 7— ray bursts observations (Milosavljevic & Nakar 
2006). Due to the very short upstream return timescale 
^ ''L/rsh|u: the particle can never stream too far ahead 
of the shock so that the turbulence is generated on small 
scales ~ 10^ — 10* cm <^ tl (Milosavljevic & Nakar 
2006). Therefore, one naturally expects in this case too 
that Fermi acceleration would be efficient, here as well, 
one needs to understand the turbulence properties before 
conclusions can be drawn on the index s. 

To summarize, we have shown in Section 2 that 
Fermi acceleration cannot operate successfully at ultra- 
relativistic shock waves if one assumes (somewhat 
naively) large-scale turbulence on both sides of the shock 
wave. The conclusions of the present discussion are thus 
more optimistic and open a wealth of new possibilities; 
in particular they suggest that the success of Fermi ac- 
celeration is intimately connected with the mechanism 
of magnetic field amplification in the shock vicinity. The 
comprehension of Fermi acceleration will eventually re- 
quire understanding the generation of the magnetic field, 
deriving the properties of the turbulence as well as char- 
acterizing the transport of accelerated particles in this 
possibly anisotropic turbulence. 

Note added: while this work was being completed, a re- 
cent preprint by Niemiec et al. (2006) appeared, report- 
ing on Fermi acceleration with small-scale turbulence. 
Although their simulations are limited to Fghiu — 10, 
these authors observe that the inclusion of small scale 
turbulence allows powerlaw spectra to emerge through 
Fermi acceleration, in good agreement with the above 
discussion. 
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